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The control of the magnetism of ultra-thin ferromagnetic layers using an electric field, rather than a current, 
has many potential technologically important applications. It is usually insisted that such control occurs via 
an electric field induced surface charge doping that modifies the magnetic anisotropy. However, it remains 
the case that a number of key experiments cannot be understood within such a scenario. Much studied is the 
spin- splitting of the conduction electrons of non-magnetic metals or semi-conductors due to the Rashba 
spin-orbit coupling. This reflects a large surface electric field. For a magnet, this same splitting is modified 
by the exchange field resulting in a large magnetic anisotropy energy via the Dzyaloshinskii-Moriya 
mechanism. This different, yet traditional, path to an electrically induced anisotropy energy can explain the 
electric field, thickness, and material dependence reported in many experiments. 

The possibility of controlling the magnetic anisotropy of thin ferromagnetic films using a static electric field E 
is of great interest since it can potentially lead to magnetic random access memory (MRAM) devices which 
require less energy than spin-torque-transfer random access memory STT-MRAM 1 " 7 . Thin magnetic films 
with a perpendicular magnetic anisotropy (PMA) are important for applications 8 ' 9 . That an interfacial internal 
electric field might be used to engineer such a PMA is also of great interest. Experiment 1011 has indeed shown that 
such a PMA might, in turn, be modified by an externally applied electric field, however the data is usually 
interpreted in terms of changes to the electronic contribution to magnetic anisotropy due to the surface doping 
induced by the applied electric field 31112 . 

The theory of the field-induced changes of the magnetic anisotropy reflecting surface doping is invariably 
developed in terms of band theory 13 " 17 . The results for both the bulk and thin films can be adequately understood 
in terms of second order perturbation theory 18 in which the relevant matrix elements of the spin-orbit interaction 
are between full and empty states. Large contributions come from regions where different d-bands (almost) cross. 
That such crossings should be close to the Fermi surface leads to a strong doping dependence in such theories. 
Nakamura et al. 19 pointed out that the strong negative applied field dependence of the PMA for an isolated mono- 
layer (ML) of Fe(001) arises directly from band splitting rather than from doping. In this case, an £ perpendicular 
to the film breaks reflection symmetry causing a large spin-orbit splitting of d-levels near the Fermi surface. As 
will be explain below, despite these important theoretical developments, a clear explanation of a number of key 
experiments is still lacking. 

Here we develop a simple analytic theory for the existence and electrical control of the PMA based upon the 
Rashba spin-orbit interaction 20 " 22 and the single band Stoner model of magnetism. We exhibit the somewhat 
delicate, but very interesting, competition between the Rashba spin-orbit fields and the exchange interaction, 
reflecting electron correlations. This theory can potentially lead to a very large magnetic anisotropy arising from 
the internal electric fields £ int which exist at, e.g., ferromagnetic/metal and ferromagnetic/oxide insulator inter- 
faces but modified by the addition of an applied electric field E ext . There is a Rashba splitting of the band structure 
leading to a quadratic, (£ int + £ ex t) 2 > contribution to the magnetic anisotropy, contrasting with a linear in £ ext 
doping effect. 



Results 

Model. This comprises a band Stoner model with the Rashba interaction added 23 : 



(1) 
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Figure 1 | (a) The electric field E = Ezis perpendicular to the ferromagnet 
surface while the order parameter direction rh, is defined by the angle 6 
relative to z. Whatever the direction of k, the Rashba magnetic field B R of 
direction kx E lies in the x — y plane, (b) The Rashba split bands of a non- 
magnetic metal. The two Fermi sheets emerge from a "Dirac point" near 
the bottom of the illustration. For the magnetic case the two Fermi sheets 
are disconnected, (c) For a perpendicular m the electron spins make a 
constant angle S to the vertical such that the projection is as in (b). The 
additional exchange splitting increases as E 2 . (d) Same but for m parallel to 
the plane. With m along the /-direction the majority and minority Fermi 
seas shift along the x-axis in opposite directions. The tilt of the spin relative 
to m is no longer a constant being zero along the x-axis and a maximum 
along the /-axis. 

where pis the electron momentum operator, 5 the order parameter, & 
the Pauli matrices and ot R = erj so E the Rashba parameter propor- 
tional to rj so , which characterises the spin- orbit coupling. The electric 
field E = Ez is taken to be perpendicular to the plane of the system, 
and ih = S/S is perpendicular to x and makes an angle 9 to the z- 
direction, as in Fig. 1(a). 

Illustrative non-magnetic example. Consider the Rashba effect in 
non-magnetic two dimensional electron gases or surface states on 
noble metals, e.g., a surface state of Au. As shown in the methods, the 
single particle energy 



H 

£ka= -—(k-ak 0 ) 
2m 

as the spin quantum number a = 

ko = maR/h 2 y and 



(2) 



:1. The momentum shift 



(3) 



identified in the inset of Fig. 1(b), reflects the single particle energy 
gain relative to zero electric field, i.e., E = 0 and hence oc R = 0. For the 



three dimensional problem e k( 



■=^-[(*- 

2m LV 



-ok 0 ) 2 + k 2 ] —E R and 



there is no equivalent momentum shift in k z . For the surface state 
of Au, E R ~ 3.5 meV 24,25 exemplifying the energy scale. 

The study of such surface states, and differences in chemical 
potentials, also helps set the scale for E. Between a metal and the 
vacuum, or dielectric, the E ~ 10 V/nm near the Cu (100) surface 26 , 



reflects the electron image potential. For the Ag/Cu(l 1 1) system 27 an 
E of similar magnitude occurs at the metal interface. The difference 
in the (1 1 1) chemical potentials 28 of Cu(4.96 eV) and Ag(4.74 eV) is 
reflected by a potential increase from Cu to Ag. The chemical poten- 
tial of Au(5.31 eV) implies a similar potential change and E between 
Cu and Au, but of the opposite sign. It is the s- electrons which 
penetrate the core and determine the spin-orbit parameter r] so . The 
example Cu/Ag/Au therefore illustrates the higher/lower potential of 
these s- electrons in the second/third transition series relative to the 
similar electrons in the 3d elements. 

There are many experiments 24 ' 25 ' 29 " 31 which put in evidence the 
Rashba splitting in two dimensional electron gases, surface states 
of noble metals, bulk layered systems, and e.g., of a surface state of 
ferromagnetic Yb. 

Magnetic case - origin of the magnetic anisotropy energy. In the 
methods it is shown the 0 dependent single particle energy, i.e., the 
equivalent of Eq. (2) for the magnetic case is: 



-,a = — [(k x - oh sin 0) 2 + k y 2 ] —E R sin 2 

- a [(j 0 s) 2 + 4 (k 2 cos 2 e + k y 2 )] 1/2 . 



(4) 



The direction of the now 9 dependent momentum shift changes sign 
as the spin index a = ± 1. These shifts also change sign with ih— ► — ih 
for a given a. This "magnetic Rashba splitting" with ih— ► — in is 
observed for the surface state of Yb 31 . 

Also in the methods the Dzyaloshinskii-Moriya (DM) and 
pseudo-dipolar (PD) contributions to the magnetic anisotropy are 
highlighted by contrasting the perpendicular and parallel orienta- 
tions of order parameter m to the plane. 

Assuming (/ 0 S) 2 > (ot R k x ) 2 and retaining the ^-dependent terms up 
to the order of E 2 in (4), we obtained our principal result: 



£an — En 



1- 



2T 



cos 2 9, 



for the magnetic anisotropy energy, with 

nr K 2 

T= — 
2m 



(5) 



(6) 



where ( ) denotes an average over the Fermi sea (see methods). The 
Rashba spin-orbit interaction produces a uni-axial anisotropy energy 
which, as in the Dzyaloshinskii-Moriya theory 32 " 35 , comprises a 
direct second order in E easy plane pseudo-dipolar interaction and 
an indirect contribution proportional to E 2 /J 0 S reflecting the com- 
petition between the first order in £, Rashba-Dzyaloshinskii-Moriya, 
and exchange fields. 

Competition between the Dzyaloshinskii-Moriya and pseudo- 
dipolar contributions. Clearly Eq. (5) implies an E 2 dependent 
PMA results if T > J 0 S/2, i.e., when the DM is larger than the PD 
term. Taken literally, the Stoner model Eq. (1), with its quadratic 
dispersion, predicts the ratio of the DM and PD contributions to the 
PMA. The result, (see methods), depends upon the spatial 
dimension. In two dimensions the PD and DM terms cancel 
although higher order terms (0(oc R 4 )) lead to a PMA while in three 
dimensions the DM term is —(4/5) E R cos 2 9 and an in-plane 
magnetisation is favoured. Lastly, a two dimensional system with a 
highly anisotropic conductivity might be modelled as a series of 
parallel one dimensional chains. For chains the DM contribution 
— (4/3) E R cos 2 9 which dominates the PD energy E R cos 2 0, 
appropriate when m is in-plane and perpendicular to the chains. 
Corresponding to the hardest axis, when m in-plane but rather 
parallel to the chains, there is neither a DM or PD contribution to 
the magnetic anisotropy energy. There is thereby a predicted electric 
field dependence of the in-plane anisotropy as seen in early 
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experiments 3 , given the large compressive strain that arises in these 
experiments. 

However, for the real problem of 3d magnets, a quadratic disper- 
sion is not at all realistic and the crystal potential V(r) must be 
accounted for, see methods. For 3d elements the wavefunction \j/ is 
well localised within the atomic sphere and the averages, e.g., 
^ 2 (fc* 2 )> and hence T, are very much increased as compared to the 
above naive estimates. In reality, the DM contribution will invariably 
lead to a PMA. 

Discussion 

The resulting anisotropy energy can be very large. The work already 
cited 24 " 27 on conducting but non-magnetic materials helps set the 
scales. The value of the scaling prefactor E R in Eq. (5) for the surface 
state of Au is ~ 3.5 meV or about 35 T in magnetic field units and 
very much larger than the typical — IT demagnetising field. If a Au 
film is polarised by contact with an ultra-thin ferromagnet, the sec- 
ond factor, 2T// 0 S, in Eq. (5) for the field inside a Au surface layer can 
be quite large ~ 5 leading to a PMA and indeed ultra-thin Fe on Au 
does have such a PMA 36,37 . Ultra-thin ferromagnetic films in contact 
with, e.g., Ru, Pd, Pt, and Ta, etc., also are found to have a 
PMA 6 ' 12 ' 38 ' 39 . 

Schematically shown in Fig. 2(a) is the potential seen by electrons 
in a free standing ultrathin ferromagnetic film. At the surface, the 
potential reflects an electron's image charge but reaches the vacuum 
level within a few atomic spacing. As already discussed, this results in 
a finite large electric field E ~ 10 V/nm at each surface but in oppos- 
ite senses. Assuming an appreciable spin-orbit coupling in the inter- 
face region, this results, in turn, in a Rashba field B R which also 
changes sign between the two surfaces for a given momentum. 
Thus, for a perfectly symmetric film, the ferromagnetically polarised 
electrons see no average field B R . This symmetry can be broken by the 
application of an external electric field as shown in Fig. 2(b). The 



electric field is increased at one surface and decreased at the other 
doubling the net effect. In contrast, for this same symmetric situ- 
ation, the surface charges are opposite and doping effects must can- 
cel. Experimentally applied fields of 1 V/nm are relatively easy to 
achieve implying a ~ 10% change in the surface anisotropy. Experi- 
ments 40 with a 1.5 nm Fe 80 B 2 o sandwiched between two MgO layers 
are perhaps closest to this situation although the thickness 1.5 nm 
and 2.5 nm of these layers are not equal. Roughly consistent with our 
estimate [see Fig. 2(c), case(i)], there is 40 an ~ 15% symmetric con- 
tribution to the magnetic anisotropy for an applied voltage of 2 V. 

Clearly the intrinsic Rashba field B R is modified when the materi- 
als adjacent to a 3d ferromagnet (F) are different. In a number of 
experiments an insulator I, often MgO, lies to one side and a normal 
metal (N), e.g., Au, Pt, Pd, Ta, or Ru, "cap" completes a tri-layer 
system. The potential, Figs. 2(d) and (e), will increase in passing from 
Fe to MgO but in passing from the F to N-layer the potential will 
either increase, Fig. 2(d), or decrease, Fig. 2(e). As discussed above, 
the relevant potential would be expected to increase, Fig. 2(d), in 
passing from a 3d element, e.g., Co, Fe, Ni, or Cu, to a 4d transition 
metal such as Ru, Pd, Ag but decrease, Fig. 2(e), for the 5d elements, 
e.g., Au, Pt, Ta. The latter case is particularly favourable since the 
intrinsic Rashba fields have the same sense and add. In addition, the 
5d elements have a larger spin -orbit coupling, resulting in a larger oc R 
and hence are more likely to produce a sizeable PMA. If the electric 
field decreases at the FI interface, the average Rashba field increases 
in the first case [Fig. 2(f)] when the effects of the surfaces tend to 
cancel and, as illustrated in Fig. 2(g), decreases in the second case 
when the inverse is true [see Fig. 2c, case (ii) and (iii)]. Experiment 12 
indeed shows an opposite field dependence for such systems with 
Pd(4d) and Pt(5d) N-layers. That the sign of the electric field con- 
tribution to the PMA reflects the N-layer whereas the field is applied 
to the opposite surface between F and I supports the current Rashba 
model. This is in stark contrast with the popular surface doping 
model 312 , for which the effects of surface doping are limited by the 



(a) 





(d) 




Figure 2 | (a) There is an electric field E in the surface region of a ferromagnet, however for a given wave vector k, the Rashba field B Ry proportional 
to k X E, has an opposite sign at the two surfaces and the average field is zero, (b) With a finite external field this symmetry is broken and there is a net 
Rashba field acting upon the electrons, (c) The gate voltage dependence of the anisotropy energy. The internal electric field causes the shift of the parabola in 
the lateral axis as indicated by V 0 for case i). For cases ii) and iii) the internal field shift is far beyond the external field range and nearly linear E-dependence 
arises, (d) The symmetry is also broken for a insulator- ferromagnet-metal sandwich. Also despite the electric field being smaller at the right surface, for a 
suitable metal, the spin-orbit coupling is larger and hence the metal interface can still dominate the net Rashba field, (e) Here the work function is 
larger for the metal than for the ferromagnet and the field for that surface is reversed. Now the Rashba fields at the two surfaces add. (f), (g) Applying positive 
gate voltages decreases the Rashba field at the insulating surface which, for this case, causes a net increase/decrease in the average Rashba field. 
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(possibly magnetically modified) Fermi-Thomas screening length. In 
reality 12 , the screening length is estimated to be much less than 1 nm, 
and much too short for there to be an appreciable doping effect of the 
Pt or Pd layers that are typically distant by a few nanometers. 

Simulating a large applied electric field £, the required asymmetry 
might be controlled in NFN tri-layers by varying in a systematic 
manner, at the mono-layer level, the thickness of one of the normal 
metal layers and by using metals with different spin-orbit couplings. 
In reality the effect of the substrate transmitted to, and through, the 
bottom normal metal will imply an asymmetry even for largish N- 
layer thickness. Indeed the PMA surface term for Au/Fe(110)/ 
Au(lll) structures does show an non-monotonic dependence on 
the top Au layer thickness 36 . Experiments 37 for Fe layers on vicinal 
Ag(OOl) and Au(OOl) surfaces and which undergo a symmetry 
breaking (5 X 20) surface reconstruction manifest an in-plane sur- 
face term reflecting this broken symmetry and which is larger for Au, 
with its stronger spin- orbit coupling, than for Ag. 

It is predicted that the surface coercivity field H c is proportional to 
(Ant + E ext ) 2 where £ int is the internal electric field corresponding to 
the zero-bias Rashba contribution to the anisotropy. Such a non- 
linear field dependence is observed, e.g., for the in-plane contribution 
for a (Ge,Mn)As/Zr0 2 surface 3 . In other experiments 10 with CoFeB/ 
MgO/CoFeB structures there is qualitative difference between the E 
dependence of the anisotropy field H c of the, "top" and "bottom", 
CoFeB layers of this three layer structure, even when they have sim- 
ilar thicknesses. The bottom layer has a larger H c and is roughly 
linear while H c becomes highly non-linear as H c — > 0 as would be 
expected as £ ext — > — £ int . 

The most direct experimental test of the model is the observation 
of the band splittings for a model Rashba system with a variable 
contact with an itinerant ferromagnet. This can result in giant mag- 
netic anisotropy (GMA) energies. For example an E R ~ 100 meV (or 
~ 1000 T) is reported in angle-resolved photoemission spectroscopy 
(ARPES) measurements 41 on bulk BiTel. For a thin film of this, or 
similar material, in contact with an itinerant ferromagnet such as Fe, 
a suitable exchange splitting / 0 S, tuned to the order of E Ry might be 
induced and a GMA will result. ARPES performed as a function of 
the direction of the magnetisation m might determine both E R and 
the momentum dependence of the exchange splitting leading to 
estimates of both the PM and DM contributions and which might 
be directly compared with magnetisation and magnetic resonance 
measurements. The electrical control of such a GMA has evident 
important application for non-volatile memory applications. There 
are clearly many more complicated embodiments of such a device. 

In conclusion, it is suggested that the Rashba magnetic field due to 
the internal electric field in the surface region of an ultra-thin ferro- 
magnet can make an important contribution to the perpendicular 
magnetic anisotropy. Such surface fields might be modified by 
application of an applied electric field. Since the internal fields at 
two surfaces tend to cancel, an asymmetry between the surfaces is 
important. Such an asymmetry is caused by different metal and 
insulator caping layers. These ideas are consistent with a large num- 
ber of experiments. 

Methods 

The non-magnetic case. This corresponds to Eq. ( 1) with / 0 = 0. It is solved by taking 
the axis of quantisation z x k to be perpendicular to the in-plane k as in Fig. 1(b). The 
eigenstates are e ik ' r \s) and H= (H 2 /2m) A; 2 —gpi B B R »<r/2, where the Rashba magnetic 
field in energy units is defined as gfi B B R = ^ a R ( — k y x + k x y) , with fi B the Bohr 
magneton andgthe g- factor, leaving the spin state |s> to be determined. There are two 
concentric Fermi surfaces. The energy splitting 2a R k = A(k/k F ), where A is the value 
for k F = (/cj?| +kF[) 1 2, with k F ^i the Fermi wave number for the spin up/down (o = 
± 1) band. For the surface state of Au, A ~ 1 10 meV while E F ~ 420 meV giving the 
E R ~ 3.5 meV cited in the text. The magnetic case. The full Eq. (1) is solved by 
defining axes such that m = S/S lies in the jy-z-plane and S = S(cos 9z + sin 9y) . The 
total field, which defines the axis of quantisation, gfi B B T = 2[(/ 0 S + a R k x sin 9)m — 
<x R kyX + ct R k x cos 9(m x x)]. It is assumed that, for a 3d ferromagnet / 0 S~ 0.5- 1.0 eV 
and gfisB R < f 0 S, i.e., the Rashba is smaller than the exchange splitting. To second 



order in gfi B B R , g[i B B T ~2(JS + a R k x sin 0)m' where/S=[(/ 0 S) 2 + 4(/c x 2 COS U + 
ky 2 )] 1 ^ 2 and where riY differs in direction from m by a small angle S where tan S ~ 
<x R (k x cos 2 9 + k y 2 ) 112 1 J Q S. The linear in k x term, a R k x sin 9, causes a shift in Fermi sea to 
give the the single particle energy Eq. (4). 

With m perpendicular to the plane, i.e., m = z (9 = 0), the exchange and Rashba 
fields are orthogonal and hence the net energy for a single electron Eq. (4) is 

tko=^k 2 -o[{hSf + {a R k) 2 ] 1 ' 2 . (7) 

The axis of quantisation is tilted by S(k) = tan - 1 — ^ away from the z-axis as shown 

in Fig. 1(c). The a = ± 1 electrons gain/lose an energy that is even in E. This arises 
from the competition of the Rashba field, perpendicular to m, with the exchange field. 
Such a competition generates a second order in E contribution to the magnetic 
anisotropy and is identified with the Dzyaloshinskii-Moriya (DM) mechanism 32 " 35 . 

Now take m parallel to the v-axis, i.e., m = y (9 = n/2). They- component of B R is 
parallel to the exchange field and is combined with the kinetic energy. The Fermi sea is 
shifted along the x-axis and lowered by E R as shown in Fig. 1(d). This energy gain 
corresponds to a pseudo- dipolar (PD) contribution to anisotropy energy 35 which 
favours an in-plane magnetisation. On the other hand, the x-component of B R , which 
is perpendicular to /oSm, gives rise to a correction to the effective exchange field. The 
direction of the moment tilts away from the y- axis in the direction perpendicular to 

the wave vector by S{k y ) =tan _1 — —f as shown in Fig. 1(d). The single particle 
JqS 

energy, Eq. (4), is now, 

h 2 

6ka = — [(k x ~ Gkof+ky 2 ] -E R 

2m (8) 

-<7[(/oS) 2 +(^) 2 ] , 

where the shift k 0 is the same as in Eq. (2) but only along the x-axis. 

The effective exchange field in Eq. (8) is smaller than that in Eq. (7) due to the 
absence of a k x term. This indicates that the overall DM contribution favours a 
perpendicular m while the PD term favours an in-plane m. This exchange field 
changes sign with a = ±1. 

Evaluation of the Dzyaloshinskii-Moriya and pseudo-dipolar contributions. Needed 

for T in Eq. (5) are the Fermi sea averages (k 2 )^ ^ = fcj^ 1 3, kp^ j / 4, and k 2 F ^ ^ j 5, 

determined analytically, for quadratic dispersion, in one, two, and three dimensions 
respectively. For an isotropic system these averages are related to J 0 S via 




which determines the ratio 2T/J 0 S in the principal result, Eq. (5), given in the text. 

Role of the crystal potential. The effects of the crystal potential V(r) are exhibited by 
considering a wave function \J/ = aK^ k+Kj ' r which is a linear combination of 
plane waves, where Kare the reciprocal lattice vectors and the a^are determined by V 
(r). While not convenient for 3d electrons, at least in principle, such an expansion in 
the true, rather than crystal, momentum states is always possible. The PD 
contribution, E R cos 2 9 is independent of the momentum k + K. However 
h 2 (k x 2 ) = \ a K\ 2 ((px + ^Q 2 ) B z' wnere ( )bz is t ne average over the first 
Brillouin zone. For 3d electrons, the average H 2 (k x 2 ), and hence T, are dominated by 
the a K for largish K. It follows T is significantly increased with the consequences 
discussed in the text. 
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